Draft version February 2, 2008 

Preprint typeset using LATpjX style cmulateapj v. 11/12/01 



in 
o 
o 

(N 
C 



> 
(N 

O 

in 
o 

^: 

Oh; 

6 ■ 



X 



HIGH RESOLUTION SPECTROSCOPY OF X-RAY QUASARS: SEARCHING FOR THE X-RAY 
ABSORPTION FROM THE WARM-HOT INTERGALACTIC MEDIUM 

Taotao Fang 1 , Claude R. Canizares 2 and Herman L. Marshall 2 

Draft version February 2, 2008 

ABSTRACT 

We present a survey of six low to moderate redshift quasars with Chandra and XMM-Newton. The 
primary goal is to search for the narrow X-ray absorption lines produced by highly ionized metals in 
the Warm-Hot Intergalactic Medium. All the X-ray spectra can be fitted by a power law with neutral 
hydrogen absorption method. The residuals that may caused by additional emission mechanisms or 
calibration uncertainties are taken account by polynomial in order to search for narrow absorption 
features. No real absorption line is detected at above 3-c level in all the spectra. We discuss the 
implications of the lack of absorption lines for our understanding of the baryon content of the universe 
and metallicity of the intergalactic medium (IGM) . We find that the non-detection of X-ray absorption 
lines indicates that the metal abundance of the IGM should be smaller than ~ 0.3 solar abundance. We 
also discuss implications of the non-detection of any local (z ~ 0) X-ray absorption associated with the 
ISM, Galactic halo or local group, such as has been seen along several other lines of sight (LOS). By 
comparing a pair of LOSs we estimate a lower limit on the hydrogen number density for the [z ~ 0) 
3C 273 absorber of n H > 4 x 10~ 3 cm" 3 . 

Subject headings: quasars: individual (PG 1407+265, PKS 2135-147, 1H 0414+009, 1ES 1028+511, 

3C 279, H 1426+268) — intergalactic medium — quasars: absorption lines — X-rays: 
galaxies — large-scale structure of universe — methods: data analysis 



1. INTRODUCTION 

Cosmological hydrodynamic simulations predict that a 
significant amount of baryons that reside in the Lya forest 
at high redshift were subsequently shock-heated to temper- 
atures between 10 5 -10 7 K, to form the so-called "Warm- 
Hot Intergalactic Medium", or WHIM (see, e.g., Hellsten, 
Gnedin, & Miralda-Escude 1998; Cen & Ostriker 1999; 
Dave et al. 2001; Fang, Bryan, & Canizares 2002; Chen, 
Weinberg, Katz, & Dave 2003; Viel et al. 2003). High tem- 
peratures ionized most of the neutral hydrogen, making it 
impossible to detect the WHIM g low redshift Lya 

forest. Thus, one of the central problems in studying the 
baryonic content of the universe is to detect these "missing 
baryons" . 

In principal, such hot intergalactic gas can be probed 
via both emission and absorption. However, the diffuse 
nature of the hot gas makes it vary hard to directly im- 
age the WHIM gas (see, e.g., Yoshikawa et al. 2003; Fang 
et al. 2005). On the contrary, much effort has been put 
into studying the absorption features produced by highly 
ionized metals in the spectra of background quasars. 

In the past few years, significant progress has been made 
in this field, thanks to the advancement in high resolu- 
tion spectroscopy provided by the new generation of space 
telescopes. A significant number of O VI absorption lines 
have been seen with the Far Ultraviolet Spectroscopic Ex- 
plorer (FUSE) and the Hubble Space Telescope (HST) 
(see, e.g., Savage, Tripp, & Lu 1998; Tripp, Savage, & 
Jenkins 2000; Tripp & Savage 2000; Simcoe, Sargent, & 
Rauch 2002). The distribution and derived properties of 
these O VI lines are consistent with predictions from sim- 
ulations. In the X-ray band, Fang et al. (2002a), Mathur, 
Weinberg, & Chen (2002), McKernan et al. (2003), and 
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Nicastro et al. (2005) reported on the detection of inter- 
vening O VII and/or O VIII absorption lines with Chandra 
and XMM-Newton. Nicastro et al. (2002), Fang, Sem- 
bach, & Canizares (2003), Rasmussen, Kahn, & Paerels 
(2003), Cagnoni et al. (2004), McKernan et al. (2004), 
and Williams et al. (2005) also reported the detection of 
z w X-ray absorption lines. These low redshift lines may 
be attributable, at least in part, to the WHIM gas in our 
Local Group. 

Starting from Chandra and XMM-Newton observation 
cycle 1, we conducted a survey of X-ray bright quasars to 
systematically study the potential X-ray absorption lines 
from the WHIM. In the first round, we observed two high 
redshift quasars (PKS 2149-306 and S5 0836+710 at z ~ 2, 
see Fang, Marshall, Bryan, & Canizares 2001) and one low 
redshift quasar (H 1821+643 at z ~ 0.3, see Fang et al. 
2002b). In this paper, we present the results from our 
second round of the survey, which includes a total of six 
Chandra and XMM targets, and a total exposure time of 
~ half- million seconds. These targets were selected largely 
based on their X-ray flux levels, which were determined by 
previous X-ray observations such as ROSAT and ASCA. 
We also selected these targets because their low Galac- 
tic absorption (except 1H 0414+009), and because their 
relatively simple spectral shape (power law) without any 
intrinsic features. Two of them were selected because their 
UV spectra showed O VI absorption lines. In principle, we 
can either take a deep look at one target, which can probe 
low column density absorbers, or observe several lines-of- 
sight to increase the total pathlength. To maximize the 
potential information we can gain, we selected the second 
approach. 

This paper is organized as follows. We present the tar- 
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gets and detailed data reduction procedures in section §2. 
We show the continuum properties in section §3. In section 
§4 we discuss how we subtract continuum and do narrow 
line analysis. Section §5 is devoted to detailed study of 
the results from narrow line analysis, and section § 6 is 
summary and discussion. 



Tabic 1: Target Parameters 



Target 




RA 


Dec 


Rcdshift 


N H 




(J2000) 


(J2000) 




(10 20 cm 2 ) 


PG 1407+265 


14 


09 23.9 


26 18 21 


0.940 


1.38 


PKS 2135-147 


21 


37 45.2 


-14 32 55 


0.200 


4.77 


1H 0414+009 


04 


16 52.4 


01 05 24 


0.287 


9.15 


1ES 1028+511 


10 


31 18.4 


50 53 36 


0.361 


1.27 


3C 279 


12 


56 11.1 


-05 47 22 


0.536 


2.21 


H 1426+428 


14 


28 32.6 


42 40 21 


0.129 


1.36 



were selected based on their strong X-ray flux obtained 
from previous ROSAT and/or ASCA observations. 



Table 2: Observation Log 


Target 


Telescope 


Obs. ID 


Obs. Date 


Duration 










(sec) 


PG 1407+265 


XMM 


0092850501 


2001-12-22 


42,062 


PKS 2135-147 


XMM 


0092850201 


2001-04-29 


59,850 


1H 0414+009 


HETGS 


2969 


2002-08-01 


52,000 






4284 


2002-08-02 


41,000 


1ES 1028+511 


HETGS 


2970 


2002-03-27 


25,000 






3472 


2002-03-28 


75,000 


3C 279 


HETGS 


2971 


2002-03-21 


107,000 


H 1426+428 


HETGS 


3568 


2003-09-08 


102,000 



2. OBSERVATION AND DATA REDUCTION 

Our sample includes six quasars with redshifts rang- 
ing from z ~ 0.13 to - 0.92. Four targets (3C 279, 
1ES 1028+511, 1H 0414+009, and H 1426+428) were ob- 
served with the High Energy Transmission Grating Spec- 
trometer (HETGS, see Canizares et al. 2005) 3 onboard 
the Chandra X-ray telescope, and two (PG 1407+265 and 
PKS 2135-147) were observed with the XMM-Newton Ob- 
servatory 4 . We list all the targets in Table 1 and their rel- 
evant observational information in Table 2. These targets 

3 See http://space.mit.edu/HETG/. 

4 See http://xmm.vilspa.esa.es/ 



2.1. Chandra Data Analysis 

The Chandra HETGS produces a zeroth order image at 
the aim-point on the focal plane detector, the ACIS-S ar- 
ray, with higher order spectra dispersed to either side (for 
ACIS-S, see Garmire et al. 2003). For all four Chandra 
observations, the telescope pointing direction was offset 
20" along +Y in order to move the zeroth order off a de- 
tector node boundary, and the Science Instrument Mod- 
ule (SIM) was moved toward the read-out row by about 
3 mm to get better ACIS energy resolution (for detailed 
instrument setups, see the Chandra Proposers' Observa- 
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Fig. 2. — Chandra MEG (red line) and HEG (dark line) spectra. Solid lines in each plot represent the best fitted models. The bottom 
panel of each plot shows the ratio between data and model. 



tory Guide, or POG 5 ). Figure la, b, c, and d display the 
zeroth-order images. We also label a size of 5" in each 
panel. 1H 0414+009, 1ES 1028+511, and H 1426+428 ap- 
pear to be point sources with core sizes less than 2 - 3 "; 
however, a large scale jet is clearly presented in the zeroth 
order image of 3C 279 (south-east direction). Detailed 
study of this X-ray jet (Marshall, Cheung, Canizares, & 
Fang 2003) will be presented in a separate paper. 

Chandra data were analyzed with the standard pipeline 
for the Chandra HETGS provided by the Chandra X- 
ray Center (CXC) 6 . We use Chandra Interactive Analy- 
sis of Observations (CIAO) V3.0. The standard screen- 
ing criteria were applied to the data. We select photon 
events with ASCA grades 0, 2, 3, 4, 6 and excluded those 
with energies above 10 keV. We also remove hot columns 
and bad pixels in each CCD chip, and take into account 
the effect of ACIS absorption. The HETGS consists of 
two different grating assemblies, the High Energy Grating 
(HEG) and the Medium Energy Grating (MEG), and pro- 
vides nearly constant spectral resolution (AA = 0.012A 
for HEG and AA = 0.023A for MEG) through the en- 
tire bandpass (HEG: 0.8-10 keV, MEG: 0.4-8 keV). The 
moderate energy resolution of the CCD detector ACIS- 
S is used to separate the overlapping orders of the dis- 
persed spectrum. We added the plus and minus sides to 
obtain the first order spectra of both grating assemblies. 
In cases where two observations were performed for one 

5 See Chandra Proposers' Observatory Guild (POG) at http://asc.h 

6 See http://asc.harvard.edu/ 

7 see http://xmm.vilspa.esa.es/ 



target (1H 0414+009 and 1ES 1028+511), we use CIAO 
tool "add_grating_spectra" to add two spectra together, 
and average the auxiliary response files (ARF). 

2.2. XMM-Newton Data Analysis 

We obtained data for PG 1407+265 and PKS 2135-147 
with all the instruments onboard XMM-Newton, and here 
we focus on EPIC (including MOS1, MOS2 and PN) and 
RGS (including RGS1 and RGS2) data. While with high 
spectral resolution RGS data will be analyzed to study 
narrow line features, we also study EPIC data to obtain 
information regarding the QSO continuum. Figure le and 
f shows the MOS images of both QSOs. 

All the EPIC observations were performed with the 
"PrimeFullWindow" and "Imaging" data mode. A 
medium filter was applied. Data were processed with the 
standard software, XMM-Newton Science Analysis Sys- 
tem (SAS) V5.4 7 . We applied the standard SAS tasks 
"emchain" and "epchain" to obtain event files for MOS 
and PN data, respectively. We selected events with pat- 
terns between and 12 for MOS and those with patterns 
between and 4 for PN. All events within central region 
with a radius of 30" were extracted for the point source 
spectra and areas between 30" and 60" for the background 
spectra. We also checked the pile-up level with the SAS 
task "epatploy" and found that there is no significant 
pile-up between 0.4 and 10 keV, so we chose this energy 

i.edu. 
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range for the subsequent spectral fitting. Finally, all the 
spectra were binned to have a minimum of 20 counts per 
bin. 

RGS 1 and 2 spectra were extracted using "rgsproc". 
This SAS task also provided a combined RMF and ARF 
file for each spectrum. To achieve the highest signal-to- 
noise, we added both RGS (including the first and second 
order spectra) together. The combined spectrum was then 
grouped channel by channel to have a constant bin size of 
0.025 A, half of the width of the RGS FWHM, for the sub- 
sequent narrow line analysis. 

3. CONTINUUM PROPERTIES OF INDIVIDUAL SOURCE 

Continuum can be fitted by a power law with absorption 
from neutral hydrogen. The flux F(E) can be expressed 
as r 

F(E) = A pl (j^j exp [-a(E)N H ] . (1) 

Here A p i is the normalization at 1 keV, T is the photon 
index, E is the photon energy, a(E) is the photoioniza- 
tion cross section, and Nh is the neutral hydrogen col- 
umn density. Table 3 lists the fitting parameters for all 
the targets and inferred fluxes and luminosities. Fitting 
was performed with XSPEC VI 1.3 8 . For each Chandra 
data set we fit with with two models: (a) a power law 
with a fixed neutral hydrogen absorption at the Galac- 
tic level, and (b) a power law with a variable hydrogen 
absorption, and we fit both HEG and MEG data simulta- 
neously (see Figure 2). The fitting is performed between 
0.9 and 8 KeV HEG, and 0.5 and 8 keV for MEG. In 
the top four targets in Table 3, the first rows show the 
fitted results from model (a) and the second rows show 
results from model (b). For the two XMM targets, since 
we cannot constrain the hydrogen column density {Nh), 
we fix them at the Galactic value and fit the data with 
model (a) only. The fitting performed with all the three 
instruments (MOS1, MOS2 and PN) simultaneously (see 
Figure 3), between 0.4 and 10 Kev. We emphasize that 

8 see http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



while from the power law fit some spectra show residu- 
als that suggest additional emission mechanisms, we are 
not trying to understand the emission processes but only 
need to fit continuum to search for absorption lines. In the 
next section, these residuals will be taken into account by 
polynomials in order to search for narrow features. Since 
in many places of the spectra we have to deal with small 
number of photons per bin, we fit the spectrum using C- 
statistics in stead of the conventional x 2 -statistics (Cash 
1976). In Table 3 and following sections, errors are quoted 
as 90% confidence unless otherwise mentioned. 
Comments on individual sources: 

• 3C 279: This source is a radio loud quasar (RLQ) 
and was identified as the first superluminal source 
(Whitney et al. 1971). Our Chandra observation, 
for the first time, revealed the jet in X-ray band 
(for detail see Marshall, Cheung, Canizares, & Fang 
2003). While both ROSAT (Comastri, Fossati, 
Ghisellini, & Molendi 1997; Sambruna 1997) and 
ASCA (Reeves & Turner 2000) observations showed 
similar photon indices (r ~ 1.8), we obtained a rel- 
atively flat spectrum with Chandra (T ~ 1.5 — 1.6). 
The fluxes and luminosities are consistent among 
these observations. 

• H 1426+428: This is a BL Lac object and 
was extensively observed with HEAO-1 (Wood 
et al. 1984), EXOSAT (Remillard et al. 1989), 
BBXRT (Madejski et al. 1992), ROSAT and ASCA 
(see, e.g., Sambruna et al. 1997), and BeppoSAX 
(Costamante et al. 2001). Our result shows that the 
fitted Njj is significantly higher than the Galactic 
value, which is consistent with results from BBXRT 
and ASCA but not with ROSAT (Sambruna et al. 
1997). The observed flux between 2 — 10 keV var- 
ied from ~ 2 to ~ 5 x 10~ n ergs cm _2 s~ 1 among 
various observations, while the photon index varies 
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Table 3: Spectral Fitting Parameters 



Target 


N H 
(10 20 cm 2 ) 


Photon Index 

(r) 


A p i a 


f b 
72-10 


£2-10 C 


C- stat/Aoi 


3C 279 


2.21 d 


1.52 + 0.01 


3.00 + 0.03 


1.4 


1.5 x 10 4t > 


1125/950 




5.0+1.0 


1.57+0.02 


3.19 + 0.07 


1.3 


1.5 x 10 46 


1104/949 


H 1426+428 


1.36 d 


1.80 + 0.01 


11.87 + 0.07 


3.6 


1.7 x 10 45 


2098/950 




9.7+0.5 


1.98 + 0.02 


14.35 + 0.20 


3.4 


1.6 x 10 45 


1456/949 


1H 0414+009 


9.15" 


2.49 + 0.02 


3.34+0.05 


0.4 


1.2 x 10 45 


1079/950 




11.88ti;f$ 


2.57+0.04 


3.57+0.13 


0.4 


1.2 x 10 45 


1068/949 


1ES 1028+511 


1.27 d 


2.11 + 0.02 


4.83 + 0.05 


0.9 


4.8 x 10 45 


1153/950 




5 - 98 -0.90 


2.24 + 0.02 


5.43 + 0.13 


0.9 


4.7 x 10 45 


1075/949 


PG 1407+265 


1.38 d 


2 24 +u - ,u 


0.44+0.01 


0.08 


4.2 x 10 45 


956/986 


PKS 2135-147 


4.77 d 


1.81 + 0.01 


1.65 + 0.01 


0.6 


6.2 x 10 44 


2400/1521 



a. Flux at 1 keV (observer's frame) in units of 10~ 3 photons cm -1 s _1 keV _1 . 

b. Flux between 2 — 10 keV (observer's frame), in units of 10 -11 ergscm~ 2 s 

c. Intrinsic luminosity between 2 — 10 keV, in units of ergss . 

d. Nh fixed at the Galactic values. 



• PKS 2135- 147: This is a typical double-lobe, ra- 
dio source (Miley & Hartsuijker 1978; Gower & 
Hutchings 1984). X-ray observations with Ein- 
stein (Wilkes & Elvis 1987) and EXOSAT (Singh, 
Rao, & Vahia 1991) reported a photon index con- 
sistent with our observation; however, ROSAT ob- 
servations (Rachen, Mannheim, & Bicrmann 1996) 
showed a much steeper spectrum, with T ~ 2.5. 
We also obtained a lower X-ray flux, compared 
with previous observations. Optical and UV spec- 
tra show strong absorption lines of O VI, N V and 
Lya at z a b s ~ z em (Hamann ct al. 1997). These 
absorption lines also lie near the center of a small 
cluster and is very close to three galaxies inside that 
cluster. It is still unclear whether these absorption 
lines are intrinsic or due to the intervening systems. 

4. NARROW LINE ANALYSIS 

Cosmological simulations predict that typical X-ray ab- 
sorption lines from the WHIM gas have line widths on the 
order of mA due to velocity broadening (see, e.g., Hcll- 
sten, Gncdin, & Miralda-Escude 1998; Fang, Bryan, & 
Canizares 2002; Chen, Weinberg, Katz, & Dave 2003). On 
the observation side, the detection of such narrow features 
relics on a careful measurement of the continuum level. 
Using 1ES 1028+511 as an example, in Figure 4 we show 
how we obtain a continuum-subtracted spectrum. 

Figure 4 shows a small portion of the raw count spec- 
trum of 1ES 1028+511 plotted against wavelength. The 
solid dark line in the top panel shows the MEG first or- 
der counts between 5 and 7 A , with a bin size of 0.02A. 
The bin size is chosen in such a way that there are at least 
two bins across the FWHM of the point spread function 
(PSF) of the instrument. In this way the bin size is 0.005A 
for Chandra HEG, O.OlA for MEG and 0.025A for XMM 
RGS 9 . The red line (labeled model 1) shows the contin- 

9 See http://asc.harvard.edu/proposer/POG/html/HETG.html and 

http://xmm.vilspa.esa.es/external/xmm_user_support/documcntation/uhb/node45.html for the width of instrument PSF of Chandra HETG 
and XMM RGS, respectively. 



slightly around 2. While BBXRT reported an X- 
ray absorption feature at <~ 0.66 keV (Madejski et 
al. 1992; Sambruna et al. 1997), we could not con- 
firm this in our Chandra MEG spectrum, consistent 
with results from recent XMM- Newton observations 
(Blustin, Page, & Branduardi-Raymont 2004). 

H 0414 +009: This target is a BL Lac object asso- 
ciated with a galaxy cluster of Abell richness of 
(Falomo, Pesce, & Treves 1993). BeppoSAX and 
ROSAT observations (Woltcr ct al. 1998) showed 
similar photon indices. However, we obtain a lower 
2 — 10 keV flux, indicating flux variation by a fac- 
tor of ~ 2. We also find a slightly high Nh, com- 
pared with the Galactic value. 

1ES 1028+511: The redshift of this BL Lac object 
has been accurately determined as z — 0.361, based 
on the measurement of two Ca II absorption lines 
(Polomskiet al. 1997). ROSAT observation showed 
a power law spectrum with T = 2.43 ± 0.20 and an 
unabsorbed flux of 6.16 x 10~ n ergs cm _2 s _1 be- 
tween 0.5 — 2.4 keV. This source was also among 
sources that listed in the ASCA Medium Sen- 
sitive Survey (Ueda et al. 2001) with a flux of 
7.77 x 10 _12 ergscm _2 s _1 between 2 — 10 keV, as- 
suming a photon index of 1.7. 



PG 1407+265 It's one of the brightest quasars at 
redshift around 1. There is weak evidence for 
Damped Lya Absorption and for Lyman Limit Sys- 
tem in its optical continuum (Lanzetta, Wolfe, & 
Turnshek 1995). A HST observation revealed in- 
trinsic narrow absorption features of O VI and Lya 
in the quasar rest frame (Ganguly et al. 2001). Pre- 
vious ASCA observation showed a spectrum with 
similar photon index (r ~ 2) but stronger 2 — 10 
keV flux (Reeves & Turner 2000). 



6 



Fang et al. 



uum from a power law plus neutral hydrogen absorption, 
with the fitting parameters adopted from Table 2.2. The 
bottom panel gives \, which effectively is the correspond- 
ing Gaussian sigma of the Poisson distribution (Gchrcls 
1986). 

While on scales larger than ~ 1 A the absorbed power 
law provides a good fit to the overall spectrum, in some 
local regions it will either over or under-estimate the ob- 
served counts due to instrumental or intrinsic fluctuations. 
Clearly, in the bottom panel of Figure 4, there are more 
red bins above zero than bins below zero, which means in 
the 5 - 7A region, model 1 underestimates the observed 
counts. The discrepancy amounts to ~ 10%. Using such 
a model could compromise the search for narrow absorp- 
tion or emission lines. To eliminate this problem, we first 
obtain the residual spectra by subtracting model 1. We 
then divide the whole residual spectra into several small 
regions, and in each region we fit the residual with a 3- 
order polynomial. To avoid significant deviation caused 
by one or two randomly high or low x bins, we ignore all 
the bins with \x\ > 3 before fitting. The selected spectral 
ranges are 2 and 30 A for MEG, 2 and 14 A for HEG, and 
5 and 35 A for RGS. The region size is 4 A, 4 A, and 5 
A, for MEG, HEG, and RGS, respectively. This method 
allows us to minimize the fluctuations at large scales, and 
still preserve features with width narrower than 1 to 2 A, 
which approximately is the region size divided by the order 
of the fitted polynomial. We find the histogram distribu- 
tion of the x obtained in this method can be well fitted by 
a Gaussian distribution. We call this method "model 2" 
and plot it in blue in Figure 4. Clearly, model 2 provides 
a significantly better fit to the observed continuum. From 
Figure 5 to 10 we showed the data and the best fitted 
results from "model 2" . 




6.0 

Wavelength (A) 

Fig. 4. — A sample of data (dark) and comparison between two 
models. Model 1 (shown in red) adopted simplily the best fitted 
model from Table 2.2, and model 2 (blue) uses an local fitting tech- 
nique. Bottom panel shows the %> the equivalent Gaussian a of 
Poisson distribution of each model. Clearly, model 2 is a better fit 
to data. 

The next step is to identify any potential absorption (or 
emission) features in the spectrum. Our search criterion 
is still based on %: we search for any feature which has 

10 see http://space.mit.edu/ASC/ISIS/ 



at least two continuous bins with \\\ > Xmin- Assuming 
the signal-to-noise ratio (SNR) in the first bin is Xl an d 
in the next bin is X2, Total SNR in these two bins follows 
X = {xi + X2)/V2- To ensure a minimum detection of at 
least 3cr, we need to have both xi an d Xi ^ — 2 for ab- 
sorption features and both xi an d Xi ^ 2 for absorption 
features. We label features identified with this require- 
ment as potential absorption or emission lines. We then 
fit the continuum-subtracted residual with a Gaussian line 
profile around the features we identified. The subtraction 
and fitting are performed with the software package ISIS 
(Interactive Spectral Interpretation System, see Houck & 
Denicola 2000) 10 

Combining all the data together, we have a total of ~ 
20,000 bins. Based on the above criteria, we find a total 
of 28 possible absorption features. Other stringent criteria 
are enforced to eliminate false detections, including that 
(1) a true feature should appear in both MEG and HEG 
spectra for Chandra targets; (2) if a feature appears in 
MEG only (at wavelength longer than 14 A), it should 
appear in both +1 and —1 orders; (3) for RGS, a true de- 
tection should appear in both RGS-1 and RGS-2; and (4) 
we should avoid regions right at the instrument features 
- these regions typically show large discrepancies between 
model and data. With these criteria no real feature is de- 
tected. We also searched for emission features and find no 
significant detections. 

Other methods have been applied to fit the contin- 
uum and detect narrow absorption features. For exam- 
ple, McKernan et al. (2004) added an inverted Gaussian 
to a physically-motivated continuum model and searched 
for absorption features by monitoring changes in fit statis- 
tics. To compare these two methods, we make a faked 
spectrum by adding narrow features to an absorbed power 
law continuum and then apply both methods. We find 
that both methods can effectively identify narrow absorp- 
tion lines with nearly the same significance. However, our 
method is more robust when the underlying continuum is 
more complex and has unknown components, and when 
there are uncertainties in continuum calibration. 

5. DISCUSSION 

5.1. Intervening IGM Absorption 

Based on the non-detection of any significant absorp- 
tion features, we can determine the minimum detectable 
ion column densities at a certain SNR level. Given the 
quasar continuum flux fx, the minimum detectable line 
equivalent width (W\) is 

W ^ = SNRx {f^T' (2) 
where A is wavelength, R = A/AA is the resolving power, 
A is the effective area and T is the exposure time. We ex- 
pect weak, unresolved features from the linear part of the 
curve-of-growth. The minimum detectable column den- 
sity (Ni) and equivalent width W\ then follow the linear 
relationship (Spitzer 1978) 

5p = 8.51 x 10~ 13 AW. (3) 
A 



Here / is the oscillator strength. 
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In Table 4 we list the minimum detectable column den- 
sities of O VII and O VIII along the line-of-sight towards 
these six targets. We adopt a SNR of 3 here. The rest 
frame wavelengths of O VII and O VIII are 21.6019 A and 
18.9689 A, respectively (Verner, Verner, & Ferland 1996). 
Since Chandra HEG has no effective area in the longer 
wavelength of both transitions, we use MEG only for the 
four Chandra targets. Since MEG has no effective area 
above 25 - 26 A, we can probe the line detection up to 
z ~ min(zi, 0.35) for O VIII and z ~ min(zi, 0.2) for 
O VII, where z% is the redshift of individual sources. 



Table 4: 3 


a Limits on 


Column Density 


Target 


O Vlll 


O Vll 


1H 0414+009 


1.87 


1.12 


1ES 1028+511 


1.44 


0.86 


H 1426+428 


1.00 


0.60 


3C 279 


1.87 


1.12 


PG 1407+265 


2.29 


1.37 


PKS 2135-147 


2.36 


1.02 



a. All column densities are in units of 10 16 cm 2 . 

Using standard cosmological models and ingredients 
from galaxy and large scale structure evolution, numerical 
simulations predict the spatial distribution of metals and 
their ionization structures (see, e.g., Cen & Ostriker 1999). 
Such simulations allow us to make quantitative predictions 
of the absorption features that could be present in the 
spectra of background sources (see, e.g., Hellsten, Gnedin, 
& Miralda-Escude 1998; Fang, Bryan, & Canizares 2002; 
Chen, Weinberg, Katz, & Dave 2003). Our non-detection 
of any significant absorption lines can, conveniently, put 
constraints on cosmological parameters and other physical 
processes adopted in the simulations. 

Rather than running complicated hydrodynamic simu- 
lations with various cosmological parameters and physical 
processes to explore the parameter space that can be con- 
strained by our observation results, we adopted a simple 
but effective analytic approach that follows Perna & Loeb 
(1998) and Fang & Canizares (2000). The basic idea is that 
all the hot gas is distributed within virialized halos that 
follow a Press-Schechter distribution (Press & Schechter 
1974). Given a certain density profile of a virialized halo, 
the probability that a random line-of-sight that passes 
through a halo of mass M with an impact parameter b 
can be calculated. Based on assumptions of metal abun- 
dance, we can then calculate the so-called "X-ray Forest 
Distribution Function" (XFDF), defined as d 2 P/(dNidz), 
the absorption line number per unit redshift (z) per col- 
umn density (iVj) for ion species i. Specifically, XFDF can 
be analytically calculated as 

d 2 P r dn d£ dl 
dNidz ~ J M dMdWidz' 1 ' 41 

where (dn/dM) is the distribution of the comoving virial- 
ized halo, E is the cross section of the halo, and I is the 
path length. 

How accurate is this analytic approach? Fang, Bryan, 
& Canizares (2002) compared both numerical and semi- 
analytic methods (with the above analytic model and a 
halo temperature profile fitted from simulations). For 
O VII and O VIII, it turns out that there is a large dis- 
crepancy in the low column density region (iVj ~ 10 12 — 



10 15 cm~ 2 ), where most lines are distributed in the fila- 
mentary structures seen in the numerical simulations but 
which cannot be described by Press-Schechter formalism. 
At the high column density end (JVf > 10 15 cm~ 2 ), the 
semi-analytic method provides a reasonablely good fit to 
results from numerical simulations. These high column 
density lines distribute typically in virialized halos with 
higher temperatures and densities than found in the fil- 
aments. Since in our study the minimum detectable col- 
umn densities are around 10 16 cm -2 , the analytic estimate 
is appropriate. 

Given the XFDF and certain cosmological parameters, 
we can calculate the expected absorption line number by 
combining observations on all the six targets. The total 
absorption line number is then 

n = £ / / aWn-Wi**. (5) 

Here Ni(J) is the minimum detectable column density of 
ion i for the jth target; Zj is its maximum redshift that can 
probed with the Chandra and XMM-Newton instruments; 
and overall summation is over all the six targets (index j). 

The parameters we plan to constrain are the cosmic mat- 
ter density il m and metal abundance Z (in units of the 
solar abundance Zq.) Changes in these two parameters 
will dramatically change gas and metal content and could 
have significantly impact on the detectability of X-ray ab- 
sorption lines. We keep all the other parameters at the 
standard values: we use the ACDM model with a dark 
energy density of Qa = 1 — fi m , the Hubble constant is 
H = 100/ikms _1 Mpc _1 with h = 0.67; the baryon den- 
sity is ilb = 0.04; and we assume the gas fraction, or the 
ratio of baryonic- to- total mass, is f gas = 0.2. 
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Fig. 11. — the expected total O VII (dark lines) and O VIII (red 
lines) line numbers with different f2 m and metal abundance Z. The 
horizontal green line shows the 90 % confidence upper limit from the 
non detection of absorption line in our sample. The vertical green 
line shows the value of Q m measured with WMAP (Spergel et al. 
2003). The green arrow indicates allowed region. 

Figure 11 shows the expected total O VII (black lines) 
and O VIII (red lines) line numbers with different fl m and 
metal abundance Z. Since we detect no absorption line, 
the horizontal green line shows the 90% confidence upper 
limit (Gehrels 1986). The vertical green line shows the 
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value of Cl m measured with WMAP (Spergel et al. 2003). 
The green arrow indicates allowed region. Clearly, we ex- 
pect to observe more absorption lines with higher values 
of metallicity and fi m . The degeneracy between metal- 
licity and Q TO can be broken by applying the much more 
accurate value of Q m measured with WMAP (Spergel et 
al. 2003): £l m /i 2 = 0.14, indicated by the vertical green 
line in Figure 11. We cannot constrain metallicity of the 
IGM based on the measurement of O VIII, i.e., the pre- 
dicted line number is still less than the 90% upper limit 
from our observation even if the metallicity is as high as 
0.3Z Q . However, O VII can put much better constraint 
on the IGM metallicity: for a WMAP value of Q m , the 
metal abundance should be smaller than ~ O.SZq. That 
O VII is more restrictive is not surprising: in the collision- 
ally ionized WHIM gas, simulations showed that O VII is 
more abundant than O VIII and should produce more ab- 
sorption lines (see Figure 13 of Fang, Bryan, & Canizares 
2002). 



5.2. Local (z w 0) Absorption 

Recently, there are a number of observations of quasars 
showing local (z w 0) X-ray absorption lines along several 
line-of-sight. These background quasars (see, Kaspi et al. 
2002; Nicastro et al. 2002; Fang, Sembach, & Canizares 
2003; Rasmussen, Kahn, & Paerels 2003; Cagnoni et al. 
2004; McKernan et al. 2004) are among the brightest ex- 
tragalactic X-ray sources in the sky (some of them are used 
as instrument calibration targets). In particular, McK- 
ernan et al. (2004) conducted a survey of nearby type 
I AGNs, and found that about half of their sample ex- 
hibit local (z w 0) absorption from H- or He-like oxy- 
gen. The high column densities of these X-ray absorbers 
(N ~ 5 x 10 15 - 2 x 10 16 cur 2 for O VII) imply the ex- 
istence of large amounts of hot gas. Given the current 
instrument resolution with Chandra and XMM-Newton, 
velocity measurements cannot distinguish its location: in 
the interstellar medium, in the distant halo, or in the Local 
Group as the intragroup medium. 

No local, z w 0, absorption lines are seen in our sample. 
We also notice none of these targets exhibit local absorp- 
tion from O VI in UV band (Sembach et al. 2003) The 
limits, from Table 4, are comparable to or below some of 
the lines detected along other lines of sight. A systematic 
study of detections and non-detections and their implica- 
tions will be presented in a forthcoming publication. Here 
we note that, whatever its location, the column density 
of local absorption fluctuates at least by a factor of two 
across the sky. More specifically, for 3C 279 our 3 a limit 
of 1.1 x 10 16 cm~ 2 for O VII is well below the detected 
value of - 1.8 x 10 16 c m - 2 for 3C 273 (Fang, Sembach, & 
Canizares 2003). These two sources are separated by an 
angular distance of ~ 10°. This corresponds to a linear 
scale of absorber structure of Al < 90(£>/500 kpc) kpc, 
for a distance to the absorber D. 11 . If the absorber is 
clumped on this scale, we can obtain a lower limit on the 



hydrogen number density of 

»// > i io '' Z(> 



D 



0.1Z C J VSOOkpc' Cm ~ 3 ' (6) 

Here we adopt a metalicity of Z<j ~ 0.1 Zq and assume 
half the oxygen is ionized to O VII. 

6. SUMMARY 

In this paper we present our Chandra and XMM- 
Newton observations of six quasars, with redshifts rang- 
ing from <~ 0.1 to <~ 0.9. Our main conclusions can be 
summarized as follows. 

1. We obtained the continuum information of these six 
targets. All of them can be fitted quite well with a 
power law with absorption from neutral hydrogen, 
although some of these targets may require addi- 
tional components to achieve a better fit. 

2. Our main purpose is to search for any narrow ab- 
sorption features in the X-ray spectra that were pro- 
duced by the WHIM gas. After applying an opti- 
mized local fitting technique to subtract continuum, 
we found a total of 30 narrow features at or above 
3ct lever. A cross check with different instruments 
excludes all features. Thus we conclude that no ab- 
sorption feature was detected in our observations. 

3. Given the size of our survey, we are able to put 
stronger constraints on cosmological parameters, 
compared with previous studies (see, e.g., Fang & 
Canizares 2000). We find the metal abundance 
of the IGM must be smaller than <~ 0.3Z Q if we 
adopted a WMAP value of £l m . 

4. We also search for local (z ~ 0) absorption lines, 
such as those that have been detected along several 
other lines of sight. No local absorption line was 
found. 

5. The limit on a (z ~ 0) absorber in 3C 279 com- 
pared with the detection of a strong absorption line 
in nearby 3C 273 indicates that the absorber has an 
angular scale of <~ 10°. If the 3C 273 absorber is in 
a cloud of this scale we estimate a lower limit on the 
hydrogen number density of nu ~ 4 x 10~ 3 cm~ 3 . 
Implications will be explored further in a forthcom- 
ing publication. 

Though several lines of evidence imply the existence of 
the WHIM gas, firm evidence for such "missing baryons" 
still remains an important challenge in cosmology. Future 
X-ray missions are very promising for detecting WHIM 
gas (see, e.g., Chen, Weinberg, Katz, & Dave 2003). Both 
Constellation-^ 12 and XEUS 13 can probe highly ionized 
metals to column densities of as low as 10 14 - 10 15 cm~ 2 . 
Such high resolution and sensitivities will reveal a true "X- 
ray Forest" in the IGM. By comparing observations with 

11 A distance of 500 kpc is likely an upper limit. It is unlikely that the 3c 273 absorber has a distance from us of > 1 Mpc since it will then be 
located even beyond our Local Group, see Fang, Sembach, & Canizares (2003) on discussion of the geometry of the Local Group and 3C 273 
absorber. 

12 See http://constellation.gsfc.nasa.gov/ 

13 See http://www.rssd. esa.int/indcx.php?projcct=XEUS. 
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numerical simulations and analytic analysis, we may even- 
tually obtain a complete theory of structure formation and 
evolution from the big bang to the present epoch. 
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